Counting Photons in the A-Experiment 
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1 Introduction 



Abstract 

Dchmclt's A-cxpcrimcnt for a three- level atom with simultaneously driven strong and 
weak transition is studied within quantum stochastic calculus approach. The statistics of 
the emitted photons is found by the method of generating functional of the corresponding 
two dimensional output counting process. In particular, the average waiting times for a 
count are calculated. 



In the eighties of last century the series of Dehmelt's papers [T][2] became the source of inspiration 
for many authors who were fascinated by the possibility of observing the fluorescence light 
emitted by a single confined atom or ion. Dehmelt suggested a very sensitive scheme for detecting 
very weak transitions in a single trapped ions (or atoms) by exploiting electron shelving effect. 
The electron shelving effect appears in atom with two transitions: one strong and one weak 
simultaneously driven. When the electron is shelved on the metastable level, the fluorescence 
light (resulting due to the intense transition) is switched off. Consequently, the fluorescence 
light emitted by a single atom exhibits periods of darkness. The length of the dark periods for 
a three-level system in the ^-configuration was obtained by Cohen- Tannoudji and Dalibard [3]. 
These results were confirmed by Barchielli [4] who gave also the physical explanation of these 
phenomena in the language of quantum stochastic calculus (QSC) Following [I] we apply 

the mathematical theory of QSC to describe the Dehmelt experiment for a three-level atom 
in the A-configuration. To find the statistics of the emitted photons we use the theory of the 
counting processes [7-11]. The generating functional approach of |11| allows us calculate the 
average waiting times for the counts. Belavkin's filtering equation for a quantum system under 
the counting observation, cf. [9] and the literature therein, enables us to study the problem in 
terms of pure posterior quantum states. 

2 The model of experiment within QSC 

Let us consider a three-level atom with two transitions: very intense |1) <->■ |0) and very weak one 
|0) <H> 1 2). We will call these transitions for convenience the "blue" and the "red" and assume 
that they are driven by two lasers. 
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Figure 1: Energy- level scheme for the A configuration. 



The hamiltonian H of the system is given by formula [12] 

(1) 



k=l 



2 ri 



2 n k (\0)(k\ + \k)(0\) + A k \k)(k\ 



where f2f~-s represent Rabi frequencies, A^s are detuning parameters. According to [3], cf. also 
[6], the dynamics of the system 'atom plus electromagnetic field' can be described by quantum 
stochastic differential equation (QSDE) if we recall to the following physical approximations: 

(i) the interaction atom-field is linear in the field operators, 

(ii) the rotating-wave approximation (RWA) is made, 

(iii) the field spectrum is flat, and the 'coupling constants' do not depend on frequency. 

Approximation (i). The interaction between an atom and electromagnetic field, in the interaction 
picture with respect to the free dynamics of the field and in the dipole approximation, takes the 
form 

— er ■ E(r, t) , 

where r is the position of the electron, and E is the quantum electric field. 



Approximation (ii). Let bj(uj), bj(uj), j G /, where I is a countable set, denote modal annihilation 

and creation operators satisfying canonical commutation relations [bj (uj) , b| (u/)] = Sji5(uj — uj'), 
where uj > is a continuous index representing energy, and j is a discrete index. The positive- 
frequency electric field operator E + (r,t) can be written in terms of bj(uj)-s as 

r+oo 

E + (r,t) = Y (27T)- 1 / 2 / F j {r,uj)e- iu]t b j (oj)duj. (2) 
- Jo 

The explicit form of the coefficients Fj(r,uj) is not important for our discussion, we stress only 
that index j also contains the direction of propagation of photons. Let ^0(^)5 V'fc(^) (k=l, 2) 
denote the wave functions of the considered states. The interaction hamiltonian in the RWA 
can be written as 
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If we introduce 



H k0 = -e\k) J Mr)r-E + (r,t)Mr)d 3 r{0\ + he. 



/f := -e / Mr)r-F 3 (r,uj)Mr)d 3 r, 



lj ■— ° J YKV J' ■*• ] 

then the interaction hamiltonian takes the form 

r+oo 

H k0 = \k) £ {^r 112 / e-^/fV)^) ^ <0| + he . 
- Jo 

Therefore, in the interaction picture with respect to free dynamics of the atom, one gets 

/+oo 
e-^/f (u + WfcJftjCc; + w fc ) dw (0| + he , 

where fc = 1, 2, denote the considered frequencies. 

Approximation (iii). We assume that fj°(co) is constant in a neighborhood of co^ and zero 
elsewhere, therefore the expression fj Q {oj + oj^) can be replaced with f^°(ojf.) (coupling constant 
independent of frequency). The range of integration can be extended from — oo to +00 according 
to the fact that the field spectrum is flat. Hence 

H k o = \k)(0\J2fj° ("k)bj(t) + he, 
3 

where 

r+oo 

bj(t) = (2TT)- 1 / 2 / e-^bjiu; + uj k ) dco . (3) 
J — 00 

Let us denote the disjoint sets of indices labeling independent field modes carrying blue and 
red photons, respectively by I±, I2. Then the coupling systematic operators are given by the 
expressions 

j ZjSi = Zj |l)(0| if j eh, , . 

3 \ z j S 2 = Zj \2)(0\ if jel 2 , w 

where Zj-s are complex coupling constants and Sj-s are lowering atomic operators. 
Let us introduce the stochastic differentials of the annihilation and creation processes: 

/■t+dt 



pt+at 

dB(t)j = B(t + dt)-B(t) = J bj(t')dt', 

rt+dt 

dB\t)j = B\t + dt)-B\t) = J b)(t')dt'. (5) 



With the help of these differentials one can represent the dynamical evolution equation for the 
unitary evolution operator for 'the system plus electromagnetic field' in the form of the Ito 
quantum stochastic differential equation (QSDE) of the form [5j [6] 

2 

dU(t) = -KU(t) dt + J2J2 ( L 3 dB )( t ) ~ rfjdBjitj) u $) > u (°) = ^ ( 6 ) 

k=ijeik 
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where 

2 



2 

k=ijei k 

and H is given by the formula (H|). With the help of one obtains 

k = iff + i(ri + zy|o><o|, (8) 

where 

A=EN 2 ' fc = l,2 (9) 
and the quantities stand for total transition rates for the blue and red transitions. 

3 The statistics of the output counting process 

Let us consider two-dimensional photon counting process 

#(«) = (io) 

with components 

t 

:= E J b}(t%(t') dt', k = l,2, (11) 

connected with blue and red photons. Here 

%{t) = U(t?bj(t)U(t), 

b){t) = U{t)\{t)U{t) (12) 

stand for the output field annihilation and creation operators [Bj, cf. also [7]. These describe 
the field after interacion with the atom within the time interval (0, t). Therefore the counting 
process is called the output counting process. We assume that all emitted photons are detected. 
If the atom is prepared initially in a pure state, then the linear Belavkin filtering equation for 
the posterior unnormalized wave function takes the form [9] 

2 

d(p{t) = -^ + Ir|o)(o|-M^)dt + ^(^-/)^(t)^(t), (is) 

3=1 

where 

r:=]TA. (14) 
k=l 

QSDE (fT5|) plays the role analogous to that of the Schrddinger equation for unobserved quantum 
system. It describes the time-evolution of a pure quantum state of the atom evolving according 
to the trajectory of the counting process. (For the case of observation process of the Wiener 
type see for instance [13] and the literature quoted therein.) 

In order to obtain the average waiting times for the counts we do not need to solve the 
equation (|13p . These can be deduced from the statistics of the counting process. It is well- 
known that the whole statistics of the counting process can be found by solving the differential 
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equation for the generating (or characteristic) functional of the process. Here we use the results 
of [11] in which the generating functional approach was applied. 

Let us consider the counting trajectory up to t, k = (j m2 ,t 2 ), • • • , (jm n ,t n ))- The 

probability density of two or more counts at the same time vanishes. It follows from that the 
detection of a photon projects the atom into one of the lower states, so the atom must be 
re-excited before the following photon is detected. To calculate the probability density (with 
respect to lTj=i dtj) of counting a photon of type j mi at time t\, a photon of type j m2 at time 
t 2 , ... a photon of type j mn at time t n , where t\ < t 2 , ■ ■ ■ , t n < t, and no other photon in the 
time interval (0, t] we use the formula |llj 

n 

p( K \t) = l[r mj (V(K\t)i>\V(K\tW), (15) 

i=i 

where 

V(n\t) = e- Kt S mn (t n ) ...S mi {h), (16) 
S m (t) = e Kt S m e- Kt , (17) 

and K is given by ([7]). 

And hence, if we assume that atom was in the state |1) at the initial moment t = 

p({j mi ,h)\t) = ^ (0 | | /c» r fe | <0|e-^ 1 1> | 2 , (18) 

where the quantity 

i£(0||fc)) = || e- K ^-^\k) || 2 (19) 

is the probability of having no counts within the time-interval (t±, t] for the atom being in the 
state \k) {k = 1, 2) at the instant t\. 

From (|18p we obtain the probability density of one count within the time- interval (to, t] condi- 
tioned by having a count of a photon of the type jo at to 

Pt((j mi , h) |(io, t )) = Pt-t ((j mi , h - t ) | \k)) k = l,2. (20) 

The form of the above expression is compatible with fact that after a count of a blue photon 
the atom is in the |1) state and after a count of a red photon the atom is in the |2) state. 
The structure of ([20]) allows us write down the probability density of detecting a blue or red 
photon at the time to + 1 and no other photons in the time-interval (to, to + 1], conditioned by 
having a count a blue photon at the instant to 

W k (t\ |1» = r fc |(0|e- K ^)|l)| 2 , k = 1,2. (21) 

If we assume that a red photon was detected at the moment to then 

W k (t\ |2» = r fc |<0|e-^*-*°)|2>| 2 , k = 1, 2, (22) 

represents the conditional probability density of detecting a blue or red photon at the instant 
to + t and no other photons in the time-interval (to, to + 1]. 
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4 Existence of dark periods 

To describe phenomena occuring in the A-experiment we will follow the approach of [3] and 
[3]. The authors of [3 J noticed that the characteristics of dark periods is controlled by the 
probabilities -P f * 0+ *(0 \\k)) for not having any count in the time-interval (to, to + t], proceeded by 
the count of photon at instant to- Due to (|20p one can put for convenience to = 0. We will also 
simplify the notation by writing P(t \\k)) instead of Pq(0 | \k)). To find these probabilities one 
should calculate probability densities Wj(t\ \k)). For this purpose let us set 

e~ Kt \k) = £>(t||fc»|j>, k = 1, 2. (23) 

3=0 

By formulas (|2T1) and (j22j) one gets 

Wj(t\\k)) = Z^M*! |fc»| 2 k = 1,2 j = l, 2. (24) 
The coefficients aj(t\ \k}) satisfy the equations of motion: 

di(t||fc» = ~n iao (t\ \k}) - iA iai (t\ \k)), 

MM \k)) = ~n iai (t\ \k)) - l -n 2 a 2 (t\ \k)) - ±r ao (t\ \k)) , (25) 

a 2 (t\ \k)) = - % -Q 2 a {t\ \k)) - iA 2 a 2 (t\ \k}) 

with the initial conditions: 

aj (0\\k)) = 6 jk , k = 1,2 j = Q, 1,2. (26) 

The solution to the system (j25[) can be obtained by using the Laplace transform technique. To 
develop formulas (f24"j) one needs only a knowledge of 

<w(*l ID) = - W, Z ' + ' A2 . ■"'+ , Z2+ ' A2 f " + , Z3+ ' A2 

2 \(z 2 -z 1 )(z 3 -zi) (zi-z 2 ){z 3 -z 2 ) (zi-z 3 )(z 2 -z 3 ) 

i„( zt+iAt _ zit | z 2 +iA 1 t zs+iAj 



a (t\ |2» = --Q 2 [ . - - e*' + e ** + ~JL ~+ e 

2 V(z 2 -zi)(z 3 -2;i) (21-^2)^3-^2) (-zi -23)^2-^3) 

where (J = 1, 2, 3) are the roots of the characteristic equation of the system (I25h . The others 
coefficients have an analogous form. 
If we assume that 

Qx > /2 2 , A » A , A » ^2 , (27) 
then the roots Zj are approximately given by 



Zl = -iA 2 + £, (28) 

2; 2 ~ 



~r - -iA x + i(r 2 - 4/2? - 4Z\ 2 - Air Ax) 1 / 2 , (29) 



z 3 ^ -ir-iiz\ 1 -i(r 2 -4/2 1 2 -4z\?-4ir^ 1 ) 1 /2 i (30 ) 

C - /2 2 {A 2 - Ax) {Ql _ ^ + 4AiA2)2 + 4rKA2 _^ )2 . (31) 
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One can observe that the real parts of the roots Zj are negative for A\ ^ A 2 . It can be proved 
that for A 1 ^ A 2 

2 p+oo 

V / Wj(t\\k))dt = 1 A: = 1,2. (32) 
j=i Jo 

This means, that the probability that at least one photon is detected in the interval (0, +oo) is 
one if A\ ^ A 2 (we have assumed that all emitted photons are detected). To prove of (J32]) let 
us notice that 

j^WMk)) = TT[(K + K^e- Kt \k)(k\e- Kh ] = -±Tr[e- Kt \k)(k\e- KH ] , (33) 



2 



Wj(t\\k))dt = 1 - t hm^Trte-^l^^le-^] . (34) 



The limit appearing on the right-hand side of (|34l) vanishes provided the real parts of the roots 
Zj are negative. Thus, (|32|) is satisfied if A\ ^ A 2 . 

Let us now consider the case of 

A X =A 2 = A. (35) 

The condition (|35|) implies that the difference in photon energies matches the energy separation 
between the two lower atomic states. One can say that this is the condition for a Raman-type 
two-photon transition between states |1) and |2). The roots Zj of characteristic equation of the 
system (|25p have now the form 



21 = -iA, (36) 

z 2 = 



-r - -iA + -(r 2 - AQ\ - AQl - AA 2 - Air A) 1 / 2 , (37) 



Z3 = 

Therefore, one can get 
and 



4 2 4 

-r - -iA - i(r 2 - aq\ - aq\ - aa 2 - AirA) 1 ' 2 . (38) 



2 ,+oo Q2 



2 



J2 i ^-(tUl))* = ^2^2' (39) 



gjj, W ^))it= W f m - (40) 

To calculate ([39]) and P0|) we have not taken into account the assumption ([2~7]h For fi\ 3> f2 2 , 
the probability that after emission of a red photon at to = at least one photon is emitted 
in the interval (0, +oo), is much less than one. This means that in this case dark periods are 
extremely long or infinite. Moreover, it is worth to notice that from (j39f) and (|40p it follows 
that there exist the possibility that the atom does not ever emit any photon. In other words, 
the population can be forever trapped. 
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Let us now analyse the situation when the condition (|35j) for equal detunings is not satisfied 
and assume that a blue photon was counted at the instant to = 0. According to (jHJ), (fT9|) . (|2"T|) 
one has 

P(0||1» = 1, (41) 

j t P(t\\l)) = (42) 
3=1 



Therefore, one gets 



Due to (1321) one can write 



p(i||i» = i - E f w 3 (t'\\i))dt> 
i=i Jo 



(43) 



2 r+oo 

p( t \\i)) = x: / 

3=1 * 



|l))dt'. (44) 



3= 

Assuming that a red photon was counted at the instant to = 0, one gets 

2 /" + 0O 

P(t||2» = £ / Wj(t'\\2))dt'. (45) 

3=1 7i 

Since the expression 1 — P(t\\k)) is the probability of at least one count in (0,t] the quantities 
p(t\\k}) := ^[l-P(t\\k))} = F short (t\\k)) + F loQg (t\\k)) (46) 

are probability densities for a waiting time for a count after emission respectively of a blue 
photon and a red one. 

Let us consider the case when a blue photon was emitted at the moment to = 0. One can write 



= 1 \?,Z v V 7T — ■ (^) 



+ iA 2 \ 2 exp(2Re(zi)i) 
4|(^2 - zi)(z 3 - Zl] 

The roots Zj are now given by the expressions (I28j) - (I31I) , We assume additionally that A\ = 0. 
From ((2ZD it follows that 

|Rezi| < |Re« 2 ,3|- (48) 

Let us set 

/■+OC 

n := / F lon Jt\\l))dt, (49) 



then by flSD and (gBJ) we obtain 

r+oo 

/ Pshort^Nl))^ = 1-71. (50) 
JO 



With the help of (ggj - 01} one gets 



W(*IH» * |r2 ^ 2| |e 22t -e^f , (51) 
Pion g (i||l)) = iI|2Re(zi)|exp[2Re(zi)t] (52) 
77 ~ (^ 1 2 -4Z\2)2 + 4r 2 Z ^2- (53) 
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Thus, from (|'27|) it follows that 77 <C 1. We introduce a time delay satisfying the following 
relation 

|2Re(z 2 , 3 « 6 « ^Re^)!- 1 . (54) 

The time-interval At between two successive counts is considered as short, if At < 6, and as 
long, if At > 9. The probability of a short waiting time after the emission of a blue photon is 

r6 r6 r+oo 

P(At<6\\l)) = / p(t\\l))dt ~ / F short (t\\l))dt ~ / F shoit (t\\l))dt = 1-/7, 
JO Jo Jo 

(55) 

and the probability of a long waiting time reads 

r+oa 

p(At>e\\i)) ~ / F long (t\\i))dt = n. (56) 

Jo 

According to (|27p . none of these probabilities depend on 6. The short waiting times are dis- 
tributed with a probability density (1 — 7I) _1 i ? s hort(i| an d the long ones with a probability 
density IJ -1 F\ ong (t\ |1)). Hence, the mean duration T snort of the short waiting intervals after the 
emission of a blue photon has the form 



' short 



z =/ tF short (t\ |1» dt ~ r/2f 2 + (57) 

The mean duration Ti ong of the long waiting intervals after the emission of a blue photon is given 
by the formula 

Hong = 77 / tF loQg (t\\l))dt = \2Rezi\~ 1 . (58) 
u Jo 

Using the assumption (|27p . one obtains in a similar way 



Q 2 T 

iWt(«||2)) ^ 



£2^ + , g .e^at 



(Zl ~ Z2)(Z3 - Z2) (Zl ~ Z 3 )(z 2 ~ Z 3 ) 

q\fw 2 

4|(z2-«i)(«3-«l)l 



(59) 



7ion g (t||2)) = ; ' ^ 2 exp[2Re(z!)t] . (60) 



It follows from (j27|) that -F s hort(^| 1 2)) is negligible. This proves that short waiting times after a 
count of a red photon are extremely infrequent. Simple calculation yields an analogous expression 
for the mean duration of the long intervals as in the case of a blue photon detected at the initial 
instant to = 0. The difference lies in the fact that in this case long time- intervals play the main 
role in describing the time evolution of the atom. 



5 Conclusion 

We have shown that there exist periods of darkness in the fluorescence light emitted by a three- 
level atom in the A configuration. But in contrast to the F-case, the atom in the A-configuration 
can stay in the coherent states for extremely long or infinite time. This statement is consistent 
with the corresponding results of |12| and [H]. The dynamics of population trapping in A-system 
has a long history. Let us stress that our main purpose was not to give new physical results, but 
to show how the theory of the counting processes can be used to find the statistical properties 
of light emitted in the A-experiment. 
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